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Two-dimensional transition metal dichalcogenide semiconductors are intriguing hosts for quantum light 
sources due to their unique opto-electronic properties. Here we report that strain gradients induced hy substrate 
patterning result in spatially and spectrally isolated quantum emitters in mono- and bi-layer WSe 2 . By correlat¬ 
ing localized excitons with localized strain-variations, we show that the quantum emitter emission energy can 
he red-tuned up to a remarkable ~ 170 meV. We probe the fine-structure, magneto-optics, and second order 
coherence of a strained emitter. These results raise the prospect to strain-engineer quantum emitter properties 
and deterministically create arrays of quantum emitters in two-dimensional semiconductors. 


Single-photon sources are crucial for a number of emerging 
quantum information processing and quantum networking ap¬ 
plication^. Of all the possible quantum emitters, solid-state 
options are leading candidates for practical applicationJ^. Un¬ 
like trapped ions or atoms, solid-state emitters are typically 
buried in a bulk, three-dimensional material which provides 
long-term stability. Further, semiconductors compatible with 
established optoelectronic processing technology enable het¬ 
erostructure devices and monolithic photonic structures for 
enhanced functionality and performance^. However, the bulk 
environment also provides obstacles and limitations for engi¬ 
neering coherent and efficient quantum light sources. First, in¬ 
teraction with phonons and fluctuating spins or c hyg es within 
the emitter’s environment can lead to dephasing33. Second, 
a high-dielectric material presents intrinsic challenges to effi¬ 
ciently extract the emitted photons into a single optical modd^. 

An intriguing alternative to host a quantum emitter is a 
two-dimensional (2D) semiconductor. Like graphene, the first 
hexagonal Brillouin zone (BZ) of monolayer (IL) transition 
metal dichalcogenide (TMD) accommodates pairs of inequiv¬ 
alent valleys. Distinctively, the IL TMDs offer energetic gaps 
at the corners of the BZ (/f-point) and strong spin-orbit cou¬ 
pling further splits the valence bands by hundreds of meV 
and the conduction bands by a few meV. This enables val¬ 
ley dependent optical selection rules which allow valley po¬ 
larization, valley coherence, and spin-valley coupling. Fur¬ 
ther, strain up to ~ 11% can be induced in these material^, 
raising the prospect for significant tuning of the electronic 
bandgajPtEl Crucially, the 2D nature eliminates the high- 
index environment which hinders photon extraction and is 
fully compatible with integrated photonics approached. 

Recently, localized excitons (OD-X) in 2D WSe 2 have 
been shown to emit non-classical light yet appear to maintain 
the general electronic and magneto-optical characteristics 
(e.g. direct band-gap, large long-range exchange interaction 
energy, large Coulomb screenii^ and large exciton g-factor) 
of the host 2D semiconductod"^. The cause of exciton 
localization has been attributed to confinement at defects 
within the electronic band-gap of the WSea. The defects 
exhibit a broad-band (~60meV) emission featurd! in IL 
and bi-layer (2L) WSe 2 that a ppear s to be energetically 
modulated near the edges of flaked^ or at wrinkled. With 
aggressive spectral filtering at these locations, the second- 
order coherence of individual emitters has been probed to 
demonstrate their quantum naturJ^StH! However, a general 


approach to obtain and potentially engineer spatially and 
spectrally isolated defects is essential for further development 
of this promising quantum photonic platform. Here we 
demonstrate that local strain gradients in the 2D crystal offer 
this capability and we take a first step towards determin¬ 
istic engineering of the emitter location and optical properties. 

Using an all-dry viscoelastic stamping procedural we 
transfer a mechanically exfoliated WSe 2 flake (from bulk, 
nanoScience Instruments product NS00182) onto a Si sub¬ 
strate pre-pattemed with etched holes with 1 pm depth and 4 
pm diameter. A confocal microscope with an objective lens 
with NA of 0.82, yielding a diffraction limited focus of ~ 460 
nm at T = 750 nm, was used for p-PL measurements with 
non-resonant CW excitation at T = 532 nm. The microscope 
used dichroic mirrors at d = 550 nm to separate the excitation 
and PL signals. The sample was placed on automated 
nanopositioners at T = 4 K in a closed-cycle cryostat with B - 
0 to 9 T superconducting magnet. All spectra were acquired 
with a 0.5 m focal length spectrometer and nitrogen-cooled 
charge-coupled device with a spectral resolution of ~40 
peV at d = 750 nm for an 1800 1/mm grating. Polarization 
dependent PL was analysed by combining a liquid-crystal 
variable retarder with a quarter wave plate and fixed linear 
polarizer. A fiber based Hanbury-Brown and Twiss interfer¬ 
ometer was used for second-order correlation measurements. 
Coincidence events from two Si avalanche photodiodes 
were recorded electronically on two synchronized input 
channels with a system timing jitter of ~ 600 ps. Efficient 
spectral filtering (resolution ~3 meV) was achieved using two 
angle-sensitive edge filters. 

An optical micrograph of the sample is shown in Fig. [T^; 
the grey region and four white disks represent the substrate 
and four etched holes, respectively, and the blue region 
corresponds to the WSe 2 flake. The two lightest contrast 
regions of the flake are assigned as IL and 2L regions. A 
portion of the IL flake is successfully placed onto a pre¬ 
patterned hole. Also observable in the contrast of the optical 
micrograph are a few wrinkles and small bumps in the flake. 
Non-resonant micro- photoluminescence (p-PL) spectra at an 
excitation power of 5 pW were obtained at the IL and 2L 
positions of the corresponding labels in Fig. [T^ as shown by 
the black and red curves in Fig. [^, respectively. As expected, 
the IL (2L) region shows a strong (weak) PL signal. The 
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FIG. 1. (a) Optical micrograph of an exfoliated flake on top of a Si substrate with etched holes (1 pm deep, 4 pm diameter), (b) Typical 
PL emission spectra for IL (black) and 2L (red) WSe 2 from the smooth and unstrained locations of the flake. The ensemble emission due to 
defect-bound excitonic complexes can be seen in both IL and 2L WSe 2 . (c-d) Color-coded spatial maps of PL with (c) integrated intensities 
in the spectral range of 680-840 nm and (d) intensities in the spectral range of 780-840 nm. (e) PL spectra corresponding to localised bright 
spots numbered in (d). The excitation power was 5 pW. The open square in (a) marks the region of spatial maps shown in (c) and (d). The two 
open circles in (c) mark the locations of the etched holes. 


IL spectrum shows two high energy peaks separated by ~ 
31meV, corresponding to the quantum-well (2D) neutral 
exciton {2D-X) and a 2D charged exciton, or trion (2D-2f^ 
The broadband emission centered at 1.68 eV f or both IL and 
2L WSe 2 is due to defect-bound excitoni^^IZU while the IL 
band-gap is direct, 2L WSe 2 exhibits an indirect bandgap 
between the conduction band minimum at S and the valence 
band maximum at K in the first BZP. The weak emission 
peak centred at 1.55 eV in the 2L spectrum is due to this 
indirect 2D-X transition. 


To investigate the spatial dependence of the PL spectrum 
on the flake, we scanned the entire region outlined by the 
yellow box in Fig. [T^. Figure [TJ; shows the spatial map of 
the integrated intensities of the emission for 680 < T < 840 
nm, where the IL region of the flake has intense PL and 
the remaining areas are dark. A one-to-one correspondence 
between Fig. and c can be seen. The two dotted circles in 
Fig. 0 mark the regions of holes underneath the flake. Note 
that the spatial map has been taken for a fixed focal plane, 
therefore the regions of the flake over the etched hole are 
significantly less bright due to alteration of the focal plane 


caused by buckling of the flake. When we generate a spatial 
map for 780 < T < 840 nm, as shown in Fig. [T|l, then the 
indirect 2T)-X emission from the 2L region of the flake is 
distinctively seen. Crucially, we also observe a few localized 
bright spots on the IL flake around the edge of the etched 
hole and a few more isolated bright spots in both the IL and 
2L regions of the flake. These spots are numbered from 1-6 
for IL and 7-9 for 2L, and their corresponding spectra are 
shown in Fig. [J. The spectra from these bright spots are 
strikingly different from the typical emission spectra of IL 
and 2L WSe 2 (see Fig.[^). They show a comb-like emission 
in the region of the defect-band and very sharp emission 
lines at much lower energies. While the comb-like emission 
has recently been reported, observation of such isolated, 
highly-detuned peaks is novel. 


To better understand the origin of the isolated highly 
red-shifted emitters, we estimate the local strain in the IL 
WSe 2 by analyzing the peak energy of the 2D-X emission at 
an excitation power of 20 [xW. The 2T)-X peak energy map for 
the IL region of the flake is shown in Fig. |^. We choose a 
reference point R on IL WSe 2 and estimate strain relative to 
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FIG. 2. (a) Color-coded spatial map of the 2D-X peak energy for IL WSe 2 which includes a few localized energy minima. The energetic shift 
of the peak with respect to reference point R was used to estimate the strain. The excitation power was 20 tiW. (h-c) (h) PL spectra 

and (c) horizontal cross-sectional A£^*^'^/strain corresponding to reference point R and energy pockets PI to P5 which correlates the effect of 
local strain gradients. The excitation power for (h) was 5 pW. 


this point. We use a calibration of -55 meV shift in the 2D-X 
peak energy for 1% uniaxial tensile strain (see Fig. S5(b) of 
Ref.G^. Note that strain estimation for the 2L is not possible 
due to a very weak indirect 2D-X emission. A strain variation 
of ~ -0.05% (compressive) to 0.3% (tensile) across the IL 
region of the flake is observed. The large portion of the IL 
region of the flake which looks flat optically (see Fig. [T^) 
shows a small strain variation (< +0.05%). More importantly, 
a few highly localized energy minima are observed. We label 
five of these positions as PI to P5. The region with the largest 
2D-X red-shift, P5, is caused by a wrinkle propagating from 
the 2L region as seen in the optical micrograph (Fig. [TJa)). 
PI corresponds to a small bump which can be seen in the 
optical micrograph. Pockets P2 - P4 are caused by bending of 
the flake around the edge of the etched hole. To characterize 
the spatial localization, we plot in Fig. ^ the change in 


2D-X peak energy for a horizontal line-cut through 

the center of the five highly strained regions. Notably, 
pockets PI - P3 show linewidths (full-width-half-maximum, 
FWFIM) approaching the diffraction limit of the confocal 
microscope, suggesting that actual strain variations are even 
more localized (< 500 nm FWFIM) than can be measured 
and are averaged out by the finite microscope resolution. 
Nevertheless, AE^^'^ (estimated strain variation) as large as 
—10 meV (~ 0.2%) is observed in the localized pockets. 

The p-PL spectra corresponding to the local minima PI to 
P5 and R are shown in Fig. |^. The spectra at the locations 
show highly red-shifted sharp emission lines together with 
the comb-like emission of the ensemble band. In the extreme 
cases (see spectrum 2, IL in Fig.[^ at location 2 in Fig.[TJl), 
a red-tuning of ~ 170 meV is observed relative to the mean 
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energy of the ensemble defect band in unstrained regions. 
Further strain-gradient correlations are found for all emitters 
in Fig. Id-e; four localized emitters occur at the edge of the 
patterned hole, the remaining occur at locations of uninten¬ 
tional strain. Conversely, the spectrum at R and elsewhere on 
the flat IL region shows only a broad band OD-X ensemble 
(not shown). This direct correlation of localized strains 
and spatial and spectral emitter isolation leads to a robust 
conclusion that local strain variations disperse the energies 
of the individual defects in the OD-X ensemble, resulting 
in comb-like emission within the nominal defect band and 
much larger red-tuning of a few defects at highly strained 
positions. Though we estimate a maximum strain of only 
~0.3%, we believe that larger strains exist at more highly 
localized spatial positions than we can measure because of 
the finite microscope resolution. 
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FIG. 3. (a) Fligh resolution PL spectrum of OD-X from location 3 
of Fig. Id. The experimental data is fit by two Lorentzian curves 
(grey line). The two peaks of the doublet are fine-structure-split OD- 
X lines and the unequal peak intensities are due to anisotropic local 
strain. Inset: Spectral fluctuations of the doublet revealed using a 
color-coded time trace of the OD-X lines. The time bin is 10 seconds, 
(b-c) The excitation power dependences of (b) integrated intensity 
and (c) linewidth of the low-energy OD-X line. The grey line is the 
saturation curve fit. The red data points in (b) and (c) mark the rep¬ 
resentative positions for the spectrum shown in (a). 


The p-PL emission spectrum at location 3 on IL WSe 2 
flake (see Fig. B) shows a sharp and spectrally isolated emit¬ 
ter at T = 782.72 nm. High-resolution spectroscopy reveals 
a doublet split by A = 726 peV with unequal intensities, as 
shown in Fig. Each component of the doublet shows 
saturation and an excitation power (Pexc) dependent linewidth 
(see Fig. and c, respectively). The solid line in Fig. 
is a saturation curve fit for a two-level system using the 
relation I = /sat (^’exc/(^’exc + ^’n)), where = 2.78 pW 
is the normalization excitation power at which integrated 
intensity of the emission peaks becomes half of the saturation 
integrated intensity (/sat)- It shows a linear power dependence 
at low /"exc and a clear saturation behavior at high Pgxc- At the 
lowest Pexc, the minimum linewidths measured are ~ 60 peV 
(FWHM). With increasing Pexc, increasing inhomogeneous 
broadening is observed. The inhomogenous broadening is 
likely caused by fluctuating charges in the environment of the 
quantum emitter at a time scale faster than the experimental 


acquisition time. Also, modest spectral wandering is observed 
at longer time scales, as shown in the inset of Fig. |^. A naive 
expectation for a quantum emitter exposed at the surface is 
severe non-radiative recombination, inhomogeneous broad¬ 
ening, and photobleaching caused by nearby surface states. 
Surprisingly, the 0D-2f states in WSe 2 , perhaps aided by the 
strong exciton binding energy, do not exhibit such deleterious 
features. Further, we And the amount of inhomogeneous 
broadening is directly linked to the non-resonant excitation 
power (Figure l^c)), suggesting that quasi or strictly resonant 
excitation could lead to minimal electric held fluctuations and 
nearly transform limited linewidths for 0D-2f states. 




FIG. 4. (a) Color-coded PL intensity as a function of polarization 
of the OD-A lines at 6 = 2.5 T in Faraday geometry. The two dot¬ 
ted vertical lines show polarization direction of low-energy (6i) and 
high-energy (6h) OD-A lines. A 40°angle between these directions 
also reflects the presence of anisotropic local strain, (b) Magnetic 
field dependence of the fine-structure split doublet, (c) Normal¬ 
ized second-order correlation function of the OD-A lines under 
non-resonant cw excitation at P = 0.5Pn- The fit (solid line) shows 
g(2)(o) = 0.17 and a lifetime of 4.14 ns. 


To investigate the nature of the doublet, we perform 
polarization-resolved and magneto-optical spectroscopy. 
Figure shows a polarization resolved intensity map of 
the doublet, revealing that the lines have unequal maximum 
intensities, are linearly polarized, and the angle between 
their polarization directions is 40°. We attribute the doublet 
to be the fine-structure of a neutral OD-A. The symmetry 
breaking of the confinement potential leads to mixing of 
the two neutral excitonic states due to the electron-hole 
spin-exchange interaction, and therefore, the emission splits 
into two lines. The energetic separation between these two 
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lines is called the fine-structure splitting Here we 

observe A = 726 [,ieV. Similar values of A were reported 
recentljD^Gll and the relatively large magnitude is likely 
due to the large long-range exchange interaction energy. 
Usually, two lines of the fine-structure split doublet have 
equal maximum intensities and are linearly polarized along 
orthogonal directions. The unequal intensities of the doublet 
lines and the non-orthogonality between their polarization 
directions of the emitter at location 3 can be ascribed to the 
presence of anisotropic strairP^. We have also performed 
polarization dependent PL for other emitters at locations 
with much less strain and find that the lines of the doublets 
show nearly equal maximum intensities and their polarization 
directions are nearly orthogonal to each other (not shown). 

The fine-structure split doublet’s behavior in an external 
magnetic field (B) is illustrated in Fig. [^. In the Faraday ge¬ 
ometry (direction of magnetic field perpendicular to the flake 

plane), A is modified according to Ab = , 

where pB is the Bohr magneton and goo-x is exciton g-factor. 
The clean fit of the data to Ab affirm that both lines of the 
doublet originate from a single neutral OD-X; the fit reveals 
gOD-x = 7.16 + 0.02. The large goo-x is not understood, but 
is in agreement with recent report xhe large goo-x 
inspires further investigations into the spin-valley degree 
of freedom in the TMD quantum emitters. Additionally, 
B also changes the central emission energy (iiAvg.) of the 
doublet due to the diamagnetic shifP^, which is given by 
Efvig. (B) - iiAvg. (0) + y (B)^, where y is the diamagnetic 
coefficient. For this OD-X, we measure y = -3.9+1.0 qeV/T^ 
(not shown). For a 3D confined system, y quantifies the com¬ 
bined contributions of confinement and Coulomb interaction 
upon application of the magnetic fieldSI. The very small y 
observed here demonstrates very strong confinement of 0D-2f. 

Finally, we observe strong anti-bunched photon emis¬ 
sion from this red-shifted emitter. Figure presents 
the second-order correlation function g*^*(T) under non¬ 
resonant cw excitation {P = Q.SPn)- The closed circles 
are measured data and the solid line is the fit using the 
relation g®(T) = 1 - where Ti is the lifetime and 

SBR-pl{\ - p) is the signal-to-background ratio (SBR). 
We obtain g®(T) = 0.17+0.02 from the fit, unambiguously 
proving quantum emission. The fit yields Ti =4.14+0.15 ns 


and the estimated SBR here is ~10. 

In summary, we have achieved resolution-limited spatial 
localization of 0D-2f with extremely low densities of ~ 1 
within a 50 nm emission bandwidth (760 < d < 810 nm) and 
tuning of the OD-X emission energy over huge range (up to « 
170 meV). We directly correlate the OD-X tuning and spatial 
and spectral isolation to microscopic pockets of large strain 
variation. These results demonstrate that strain engineering is 
a viable approach to obtain spatially and spectrally isolated 
quantum emitters in 2-D semiconductors. The passive exciton 
emission tuning observed with the 2D quantum emitters 
is significantly larger than the state- of-the-art for quantum 

emitters in bulk semiconductorJ^llzejgl] However, there is 

significantly more potential for strain engineering OD-X 
states. First of all, massive strain gradients are possible in 
2D materialJ^ and a number of approaches for both static or 
in-situ tunable uniaxial and biaxial strain engineering can be 
considered^. Secondly, the affect of strain on the electronic 
and optical properties of OD-X states beyond their emission 
energy has yet to be investigated. For instance, in-situ strain 
tuning can enable engineering of the carrier confinement 
potenti aP^, the permanent dipole moment, the fine-structure 
splittinJ^Hmil and the spin propertied^ of localised exciton 
states. Notably, such quantum emitt ers c ould underpin a 
hybrid semiconductor-atomic interfacd^^l^Zl Finally, rather 
than suspending the flake over a large hole, by placing the 
flake over a nanostructured surface such as a periodic array 
of steps or holes, one can envisage periodic arrays of strained 
quantum emitters. 
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